Determination of a,^ from Gross-Llewellyn Smith sum rule by 
accounting for infrared renormalonB 



Carlos Contreras^, Gorazd Cvetic^, Kwang Sik Jeong^ and Taekoon Lee^ 
^Dept. of Physics, Universidad Tecnica Federico Santa Maria, Valparaiso, Chile 
"^Dept. of Physics, Korea Advanced Institute of Science and Technology, Daejon, Korea 



Abstract 

We recapitulate the method which resums the truncated perturbation series 
of a physical observable in a way which takes into account the structure of 
the leading infrared renormalon. We apply the method to the Gross-Llewellyn 
Smith (GLS) sum rule. By confronting the obtained result with the experi- 
mentally extracted GLS value, we determine the value of the QCD coupling 
parameter which turns out to agree with the present world average. 



The infrared (IR) renormalons are singularities at positive values of the Borel variable 
b in the Borel transforms of QCD physical observables. These renormalons represent the 
contributions of loop integrations in the low energy regimes [0]. They represent obstacles for 
and lead to ambiguities in the Borel integration. This integration then gives discontinuities 
(cuts) of the observables in the complex coupling parameter plane. It is thus important to 
eliminate these unphysical cuts from the Borel integration. This approach was proposed in 
and applied to the GLS sum rule in 0. Here we recapitulate the method and the results 
of and emphasize some physical aspects. The GLS sum rule is the quantity 



M3(Q') = 3 (l - A(g2)) = £ dx Fs{x; 



l + O 



(1) 



where = {F^^ + F^^)/2 is the nonsinglet structure function in the uN DIS, niN is the 
nucleon mass. Theoretical knowledge about the GLS A((5^) is the following. The truncated 
perturbation series (TPS) up to [^(Q) = c^siQ)/'^] is known [Q 

A(Q2) ^ a(l + ^^a + W2a^) , (2) 

as is also the structure of the first IR renormalon singularity of its Borel transform 

BT(6) = ^^^y^[l + «:i(l -&) + ■■■] . (3) 



*Partly based on invited talk given by G.C. in WG3 of NuFact02, July 1-6, 2002, Imperial College, London 
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Here v = (/5i//5o — 72)//5o, where j3o and j3i are the one- and two-loop QCD beta function 
coefficients and 72 is the one-loop coefficient of the anomalous dimension of the twist-four 
{d = 2) operator appearing in the Operator Product Expansion (OPE) for A((5^). When the 
number of active quark flavors is ?2/ = 3, we have /?o = 9/4, /?i = 4, and 72 = 8/9 (the latter 
value obtained in P), giving v = 32/81. The TPS of BT(6) in powers of b is known to 
as a consequence of (^. A consequence of the IR singularity (^ at 6 > 1 is that the Borel 
integral Abi{z), which for the complex coupling z = \z\ exp(i0) is defined as 



+00 exp(i</ 



db exp 



BT(6) 



(6 = \b\e^'^) , 



(4) 



has a discontinuity (cut) at the positive real axis z = a{Q) > 0. Namely, for z near the 
positive real axis, it can be shown from (ffl) that 



BT(6) . 



/?oa(Q). 

This quantity is not even real; the discontinuity shows up in its imaginary part 



(5) 



Im. A-Qi{a{Q):Lie) = -5-Ini / cifeexp 



c 



(1-6) 



^-^ L-l/(/3oa(Q)) 



a(Q) "V{—v) sin(7rz/) [1 + KiPoau + ■ ■ 



(6) 



Because {—a^ie) '^ = a '^[cos(7ri/)±i sin(7rz/)], (|^) leads to 



Abi(z) 



^l/(/3oa{Q)) 



+ ■ 



+ Abi(^), (7) 



where A-bi{z) is a function without cuts in the complex 2;-plane since the ffist expression on 
the right-hand side absorbs the renormalon-induced cut at z > 0. This expression, when 
z = a{Q) is real, has the Q-dependence Q~'^as{Q)'^'^^^°[\- + 0{as)\, the same as the twist-four 
term in the OPE for the GLS sum rule 0. The Im part of this expression, i.e. expression (^, 
can be identified as the Im contribution from the leading IR renormalon (^. The premise 
of the method is that the full ffist expression on the right-hand side of (|^), i.e., including its 
real part, represents the unphysical cut contribution which is to be eliminated. This leads 
to the final result for the resummed GLS value 



A{a{Q)) = ABi(a(Q)) = [Re =F cot(7rz/)Im] / o?6exp 

■J ±ie 



Poa{Q) 



R{b) 



[i-bY+''' 



(8) 



In the Borel integration here, the exactly known IR renormalon singularity has been factored 
out explicitly. Function R{b) = (1 — bY^^BTip), whose TPS up to is known exactly due 
to (0), has much weaker singularity at 6 = 1 than BT(6). The function R{b) has, in principle, 
singularities from other renormalons, at 6 = —1, ±2, ... To accelerate the convergence in the 
TPS for R further, we introduce in integral (^ the change of variables b \—>- w according to 
the conformal mapping 



wib) = (^VTTb - ^1 - 6/2) (^ViTb + ^1 - 6/2) 



-1 



(9) 
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which maps all these renormalon locations onto the unit circle \w\ = 1, except for the leading 
IR renormalon w{b = l) = 1/3 whose structure has been factored out. Similar applications of 
conformal mappings have been applied in other methods in Refs. §,0]. Further, expansion 
of R{b{w)) in powers of w suggests that the N^LO coefficient (at a^) in the GLS TPS (@), in 
MS and with renormalization scale fi = Q, is W3 = 158 ± 30, somewhat higher than the TPS 
PMS prediction ws ^ 130 g. Using u;3 = 158±30, we use for R{b{w)) in (D the TPS up to 
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Expression (|]) represents the massless part of the GLS sum rule. For ^ 2-3 GeV^, 
where reasonable experimental data are available, the contributions of the massive c-quarks 
0] represent corrections contributing about 3-4% to A((5^). Nuclear effects of the iron 
target ^ contribute about 1% to A{Q^) [Afc(Q^) ^ 4 • lO-^GeVVQ^]- The "light-by- 
light" contribution in the GLS TPS (Q) (it is ~ a^) is kept apart, it contributes to A(Q^) 
less than 0.5%. 

We work in the MS renormalization scheme. The renormalization scale fi was fixed 
according to the principle of minimal sensitivity (PMS), obtaining /x^ ^ 3.3Q^. Experimen- 
tally extracted values of GLS A((5^) at = 3.16 GeV^ are given by the Fermilab CCFR 



Collaboration [11 



A(g2 = 3.16GeV2) = 0.149 ±0.039 , (10) 

where the nuclear corrections were subtracted out. When we confront (|10|) with the theo- 
retical GLS expression [Eq. (|^) with the conformal mapping (P); plus c-quark corrections; 
plus "light-by-light;"], the obtained QCD in MS is 

a,{Q) = 0.305l[!:J^^(exp) ± 0.006(th) , as(Mz) = 0.1167^|]:[^}?^(exp) ± 0.0008(th). (11) 

Using the central CCFR data value of A((5^) at Q^ = 2.00 GeV^, we obtain almost the same 
central value: as{Mz) = 0.1166. The values of as{Mz) obtained are in agreement with the 
present world average. The experimental uncertainties are large. The theoretical uncer- 
tainties are dominated by the ambiguity with respect to the change of the renormalization 
scheme. 

The term proportional to cot(7rz/) in (^, which eliminates the unphysical cut, was im- 
portant numerically as it decreases the value of A((5^) by about 10%, thus increasing the 
extracted value of as{Mz) by about 0.003. The CCFR Collaboration by evaluating the 
NNLO TPS of A(Q2)^ corrected by a term ocl/Q^ as estimated by QCD sum rules, obtained 
from the central values of the data as{Mz) ~ 0.113-0.114. An earlier similar (TPS-PMS) 
evaluation @ gave as{Mz) !^ 0.115. 
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